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1 Introduction 



The purpose of this paper is a further study of the Kazakov-Migdal lattice gauge theory []J 
which is defined by the partition function 

Zkm = / [[dU^x) [[d<$>(x)e V * J. (l.i) 

Here the field $(se) takes values in the adjoint representation of the gauge group SU(N) 
and the link variable U^x) is an element of the group. The lowest coefficient of the 
expansion of the potential V[<&\ in the lattice spacing, a, 

V[$] = m^ 2 + ... (1.2) 

is identified with a bare mass parameter while the others play the role of couplings of 
self- interaction of the field $. 



The recent interest in the model (LI) is due to the following reasons. As it is proposed 
in Ref. the model (1.1) induces QCD in the continuum limit which should be obtained, 
as usually in lattice gauge theories, in the vicinity of a second order phase transition. 
This limit should be reached by approaching the couplings of the potential ( |1.2j ) to the 
corresponding critical values. As was proposed in Ref. PJ this critical point should be 
identified with the one separating the strong and weak coupling phases of the model ( |1 . 1| ) . 
On the other hand the model (|1 . 1|) is potentially solvable in the large- N limit with different 
distributions of eigenvalues of the matrix $ in the two phases [|l], [2[ . To obtain a confining 
continuum limit, one should approach the critical point from the site of strong coupling 
phase while approaching from the weak coupling phase seems to result in deconfining 
Higgs phase. 

However, soon after this scenario of inducing QCD had been proposed, it was pointed 
out by Kogan, Semenoff and Weiss H that the model ( |1.1| ) possesses an extra local Z N 
symmetry which leads in the strong coupling phase to local confinement when quarks can 
not propagate even inside hadrons. To escape this physically unacceptable picture, it was 
assumed this symmetry to be spontaneously broken at the point of the strong to weak 
coupling phase transition. A similar scenario has been proposed by Khokhlachev and the 
author when normal confinement restores after the large- N phase transition with occurs 
before the one associated with continuum limit {i.e. within the strong coupling phase in 
our terminology). This conjectured phase looks similar to the weak coupling phase of 
standard lattice gauge theories. To avoid terminological confusions I shall refer it as the 
intermediate coupling phase. 

The problem of whether the large- N phase transition occurs according to this scenario 
is a dynamical one and can be studied using standard methods of lattice gauge theories. 
In particularly, the mean field method has been applied in Ref. [Q. The result for the 
case of the quadratic potential V[<&] in negative — no first order phase transition occurs 
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in the strong coupling phase. This conclusion coincides with the one made on the basis 
of the exact solution found for the quadratic potential by Gross JJ. 

While the above scenario failed for the quadratic potential for this reason, a possibility 
of an alternative 'stringy' continuum limit has been conjectured for this case by Kogan, 
Morozov, Semenoff and Weiss when the lattice spacing is taken to be iV-dependent 
and should approach zero simultaneously with N — > oo in a special way. Such a 'stringy' 
large- N limit differs from the 't Hooft topological expansion of QCD which is dictated by 
the known dependence of the coupling constant on N prescribed by asymptotic freedom 
at small distances. For this reason, all the standard large-iV technology like factoriza- 
tion, saddle point equations, etc. are not applicable to the 'stringy' large-iV limit. In 
particularly, the exact N = oo solutions of Refs. [0, || can not be applied as well. 

The latter exact solutions have been obtained by solving the 'master field equation' 
derived by Migdal under the assumption that the path integral over is saturated 
as N — ► oo by a single x-independent saddle point configuration $ s — the master field. 
While it was pointed out |j| that such a master field does not contradict to all our 
knowledge about the large- N limit, a mechanism of its appearance was mysterious. 

In the present paper I consider this problem from the viewpoint of the large- N reduc- 
tion which had been first advocated by Eguchi and Kawai |7J for lattice gauge theories at 
N = oo. The large- N reduction states that the model on an infinite lattice is equivalent 
that at one point (a plaquette in the case of lattice gauge theory) so that the space-time 
degrees of freedom are eaten by the internal symmetry group. My idea would be to iden- 
tify the master field $ s with a the saddle point configuration of a one-matrix model which 
appears after the reduction. I consider both quenching || and twisting || prescriptions 
of the large- N reduction and argue in Section 2 that while they correctly reproduce the 
perturbative expansion of the Kazakov-Migdal model, any dependence on the quenched 
momenta can be absorbed by a (nonperturbative) gauge transformation. 

To justify this reduction procedure, I derive in Section 3 loop equations for the 
Kazakov-Migdal model. I show that in addition to the adjoint Wilson loop 

W A (C) = ^(\trU(C)\ 2 -l)) (1.3) 

where the average is understood with the same measure as in Eq. ( [Lip , the objects of a 
new kind 

G(C xy ) = (£tr mx)U(C xy My)U(C yx )) ) (1.4) 

emerges in the loop equations. This has a meaning of the average for an open loop which 
is made gauge invariant by attaching scalar fields at the ends. For the quadratic potential 
I obtain in the large- N limit the closed set of two equations. 

In Section 4, 1 discuss properties of the loop averages both in reduced and in unreduced 
cases and show explicitly of how they coincide to the leading order of the large mass 
expansion. I speculate as well on the properties of the intermediate phase and argue that 
it should resemble the standard Wilson lattice gauge theory. 
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An exact solution of the loop equations in the strong coupling phase is found in 
Section 4 for the quadratic potential at any number of dimensions D. An explicit formula 
for the average ( |1.4j ) is given by Eq. ( |5.7| ) below. While this solution agrees with that by 
Gross I do not make any assumptions about the master field to find it. Moreover it 
is a first example of exact calculations of extended objects in Kazakov-Migdal model. I 
extend the solution to a more general potential 

NtrV[$] = m 2 iVtr$ 2 + iV 2 ,/(i tr $ 2 ), (1.5) 

where / is an arbitrary function, of the type studied recently in the matrix models [fLOf 
and show that no large- iV phase transition occurs in the strong coupling phase for this 
potential as well. 



2 Large- N reduction 
2.1 Scalar field 

The idea of large- N reduction was putted forward by Eguchi and Kawai |7j who showed 
the Wilson lattice gauge theory on a D-dimensional hypercubic lattice to be equivalent at 
iV = oo to the one on a hypercube with periodic boundary conditions. This construction 
is based on an extra (Zjy) I) -symmetry which the latter theory possesses to each order 
of the strong coupling expansion but is broken in the weak coupling region. To cure 
the construction at weak coupling, the quenching prescription was proposed by Bhanot, 
Heller and Neuberger || and elaborated by many authors (for a review, see ||11|1). An 
elegant alternative reduction procedure based on twisting prescription was advocated by 
Gonzalez- Arroyo and Okawa ||. An extension of the quenched Eguchi-Kawai model 
to the case of hermitian matrices was proposed by Parisi [12| end elaborated by Gross 
and Kitazawa [13j while that of the twisting prescription was advocated by Eguchi and 
Nakayama [|14j] and has been discussed recently by Alvarez-Gaume and Barbon [[10] in the 
context of D > 1 strings. Let me first briefly review these results which allow to reduce 
the partition function of self-interacting matrix scalar field on the infinite Z)-dimensional 
lattice at iV = oo to a hermitian one-matrix model in an external field. 

For a pure scalar theory whose partition function is defined by the path integral similar 
to ( |1.1|) but without gauging: 

f -i-r , s y Ntr ( -VWx)]+Y *(i)*(a;+u) ) 
Z = n d$ ( X ) e ^ ' ' ( 21 ) 

J x 

the quenched momentum prescription is formulated as follows. One substitutes 

$(x) -> S{x)$S\x) (2.2) 

where 

[S(x)] tl = (2.3) 
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is a unitary matrix which eats the coordinate dependence. The averaging of a functional 
with the same weight as in Eq. ( |2.1| ) can be calculated at iV = oo by 

D N dk? 

n=\ i=i 

where the average on the r.h.s. is calculated for the quenched reduced model defined by 



w*)i> - L n n (nsM^wi)^ ^ 



(2.5) 



the partition function |12 

z QRM = J d * e -»*™ + »^K D -^™K-W 

which can be obtained from the one (|2.1|) by the substitution (|2.2| ). 

Since N — > oo it is not necessary to integrate over the quenched momenta in Eq. 
The integral should be recovered if fcf s would be uniformly distributed in a D-dimensional 
hypercube. Moreover, a similar property holds for the matrix integral over $ as well which 
can be substituted by its value at the saddle point configuration $ s : 

-> F[S(x)$ s S\x)} . (2.6) 



This saddle point configuration was referred as the master field [15 



An alternative reduction procedure is based on the twisting prescription. One performs 
again the unitary transformation (|2.2j ) with the matrices S(x) being expressed via a set 
of D (unitary) N x N matrices by the path-dependent factors 

S(x) = P 1] r M . (2.7) 

i&c x 

The path-ordered product in this formula runs over all links I = (z, fi) forming a path 
C xoo from infinity to the point x. The matrices are explicitly constructed in Ref. 
and commute by 

Y fJL T v = Z^TvTfi (2.8) 
with Z^y = Z\ being elements of Z^. 

Due to Eq. ( |2.8| ), changing the form of the path multiplies S(x) by the abelian factor 

Z(C)= 1] ^{U) (2.9) 

neS:dS=C 

where (fi, v) is the orientation of the plaquette □. The product runs over any surface 
spanned by the closed loop C which is obtained by passing the original path forward and 
the new path backward. Due to the Bianchi identity 

II Z, u (a) = i (2.10) 

where the product goes over six plaquettes forming a 3-dimensional cube on the lattice, 
the product on the r.h.s. of Eq. (|2.9| ) does not depend on the form of the surface S and 
is a functional of the loop C. 
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It is easy to see now that under this change of the path one gets 

[S(xMS\x)] kl - \Z(C)\ 2 [S(x)} tJ [S^x)] kl (2.11) 



JU I 

,2 



and the path-dependence is cancelled because \Z{C)Y = 1. This is a general property 
which holds for the twisting reduction prescription of any even (i.e. invariant under the 
center Z N ) representation of SU(N). 

For definitiveness one can choose 

s{x) = r^r^r^q 4 (2.12) 

where the coordinates of the (lattice) vector 3j til" 6 measured in the lattice units. 

For the twisting reduction prescription, Eq. (|2.4j) is valid providing the average on 



the r.h.s. is calculated for the twisted reduced model which is defined by the partition 



function [14 

Z TRM = Jd$e- NtrV[ * ]+N ^ tTV »* r > . (2.13) 

Now the perturbation theory for the unreduced model fl2.1|) is recovered due to an explicit 
dependence of T M on momenta. 

2.2 Kazakov— Migdal model 

Analogous quenching and twisting reduction prescriptions hold for the Kazakov-Migdal 
model as well. To make consideration similar to that of the previous section, let us 
introduce for the Itzykson-Zuber-Metha integral the notation 

![$,*] = / dUe"**™* = ^l^Ti ( 2 - 14 ) 



A[$]A[tt] 

where dll is the Haar measure on SU(N) while $j and stand for eigenvalues of the 
matrices $ and respectively, with A[$] = Tii<j(®i ~ $j) being the Vandermonde 
determinant. This formula implies |IJ the following representation of the partition func- 
tion ( |1 . 1|) in terms of the path integral over $(x): 

Zkm = f LI d^{x) e " " tr v[ * ix)] J] , $(x + n)] . (2. 15) 

X X,[l 

It is instructive to refer Eq. (|2.15| ) as the matrix model representation of the partition 
function fll.lj) . 

One can apply now to Eq. ( [2.15| ) at = oo the reduction prescription described in 



the previous section which results in Eq. (2.4) (or ( |2.6|) ) with the reduced models defined 



by the formulas like Q2.5p for the quenching prescription and or like ( |2.13| ) for the twisting 
prescription. However, since the Itzykson-Zuber-Metha integral depends only on the 
eigenvalues of and $(x + //), it does not depend actually on S(x) and S(x + fi): 

I[S(x)$S*(x), S(x + /j,)$S^(x + /x)] = /[$, $] . (2.16) 
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Therefore, the dependence on the quenched momenta or on the matrix which is con- 
structed from momenta is cancelled. 

The averages in the Kazakov-Migdal model can now be calculated according to 
Eq. ( p.6| ) via the master field. For the free energy itself, one gets 

^j\ogZ KM = 21og(A[$ s ]) -NtrV[$ s ] + D log/[$ s , $ s ] (2.17) 

where Vol. stands for the number of sites on the lattice. 

Some comments concerning the proposed reduction procedure of the Kazakov-Migdal 
model are now in order: 

• While the reduced model involves no explicit dependence on the momenta entering 
either S(x) or T M , this dependence survives when calculating averages of x-dependent 
quantities according to Eq. ( |2.6| ). The point is that the invariance which remains 
after the reduction is solely global SU(N): $ s — > Q&sfl^, which can cancel the 
matrix S(x) at only one point (this corresponds to the translation invariance of 
averages in the unreduced model (|L1|)). 

• Eq. Q2.16| ) is a consequence from the invariance of the Haar measure in Eq. Q2.14T) 
under multiplication by a unitary matrix from the left and from the right separately, 
i.e. S(x) and S(x + /x) which appear after the substitution (|2.2|) can be absorbed 
by a gauge transformation. While only SU(N) gauge transformations are allowed, 
this is enough for our purposes since the U(l) part is again cancelled in the bilinear 
expressions of the form [S(x)]ij[S^ (x)]ki- 

It is instructive to see of how the planar graphs of the original model ( |1 . 1| ) are repro- 
duced by the proposed reduction procedure employing the known results for the reduction 
of scalar field reported in the previous section. To this aid, let us calculate the partition 
function (|1 . 1| ) first integrating over scalar field: 



Zkm= / II dU M e ~ Sind[UKx)] , (2.18) 

where the induced action for the gauge field U^x) is defined by the integral over $(x) in 
Eq. (0): 



g -SindiU^x)} 



I shall refer Eqs. Q2.18j ) and ( 2.19 ) as the gauge field representation of the Kazakov-Migdal 
model. 



The required planar graphs of the model fll.l]) can then be obtained in two steps. The 
first one is to calculate the induced action employing the reduction prescription to the 
integral over scalar field on the r.h.s. of Eq. ([2.19| ). The second step would be to average 
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over the gauge field according to Eq. (|2.18|) . Since I do not apply the reduction procedure 



for the gauge field, it would be enough to show that the two induced actions coincides 
perturbatively. 



For the r.h.s. of Eq. ( |2.19| ), one can apply the formulas of the previous section which 



gives for the quenched momentum prescription: 

SiUU^x)] = NJ2 tr ($ 8 St(x)Ur(x)S(x + ii)Q a St(x + M(x)S(xj) . (2.20) 



X,fl 



Expanding the matrices U^{x) around the unity and using standard rules [11| to obtain 
diagrams for the scalar field, one recovers the correct diagrammatic expansion for the 
induced action Si n d[U^(x)}. 

Notice that when one integrates over U^{x) as discussed above, these diagrams are 
absorbed into gauge degrees of freedom. This looks precisely like the mechanism of recov- 
ering planar graphs by the master field which has been proposed in Ref. JIJ on the basis 
on the conjecture about the translationally invariant master field. I have related therefore 
this scenario to the large- N reduction phenomenon and the master field of Ref. to the 
saddle point matrix $ s . 

Up to now I did not discuss what equation determines $ s . This equation can be derived 
by averaging the (quantum) equation of motion for $ field in the Kazakov-Migdal model, 
written in the form (|2.15|) , according to the prescription ( |2.6| ). Introducing the spectral 
density p s (X) which describes the distribution of eigenvalues of the matrix <3> s and using 
the remaining global gauge invariance, one gets 



dx _ Ps (x) dV(X) D d 5logI[p,p s 



X-x dX N 2 dX Sp{X) 



[2.21) 



P=Ps 



which coincides with the saddle point equation for the coordinate-independent master 
field of Ref. 0. 

I did not presented a direct proof of the reduction prescription (|2.6| ) with $ s given 
by Eq. (|2.21|) . The above perturbative arguments work only in the weak coupling region 
where the perturbative expansion is relevant and are not applicable in the strong coupling 
region which is separated by a phase transition. I hope that loop equations which are 
derived in the next section could be useful for this purpose. 



3 Loop equations 

As is pointed out in Ref. |4j], the simplest gauge invariant objects in the Kazakov-Migdal 
model are the adjoint Wilson loops which are defined by Eq. ( |1.3|) where the average is 
understood with the same measure as in Eq. The nonabelian phase factor U(C XX ) 

which is associated with a parallel transport from a point x along a closed loop C xx is 
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defined by the path ordered product 



U{C XX ) = P J] U^z) (3.1) 

where I stands for the link (z, /i). 

One can derive the set of loop equations satisfied by these quantities quite similarly to 



the case of standard lattice gauge theory (for a review, see |T6|]) performing an infinitesimal 
shift 

U fl (x)^(l + ie lt (x)U lt (x)) (3.2) 

of the link variable U^x) at the link (x, fi) with e being an infinitesimal hermitian matrix 
which leaves the Haar measure invariant. Since the plaquette term is absent in the action 
of the Kazakov-Migdal model, it does not appear on the l.h.s. of the loop equation. 
However, a new term associated with the interaction between gauge and scalar fields 
arises on the l.h.s.. The r.h.s. of the loop equation satisfied by Wa(C) looks at large N 
very similar to the standard one. The equation reads schematically 

(^tr (§{x)U{C xx )U,{x)${x + ii)Ul{x)U\C xx ] 

(itr (®{x)U\C xx )U^x)®{x + v)Ul{x)U{C xx ) 

E r,(l)5 xz W A (C xz )W A (C zx ) (3.3) 

where r M (Z) stands for a unit vector in the direction of the link /. 

A conceptual difference between Eq. ( |3.3| ) and the standard loop equation of manycolor 
QCD is that the former one is not closed. While the r.h.s. of Eq. ( ^.3| ) involves the same 
quantity Wa, a new gauge invariant object of a generic type ( |1.4|) emerges on the l.h.s.. 
One should derive, therefore, an equation satisfied by this quantity. 

The corresponding equation results from the invariance of the measure over $ under 
an infinitesimal shift 

-> (3.4) 

of at the given site x with £(x) being an infinitesimal hermitian matrix. While to 
close the set of equations one has to consider averages of the type (|1.4|) with arbitrary 
powers of and $(y), drastic simplifications occur for the quadratic potential when 
the resulting equation reads 

D 

2m 2 G(C xy ) - E G(C (x+ , )x C xy ) = 5 xy W A (C xy ) . (3.5) 

fi=-D 

Here the path C( x+fJ ) x C xx is obtained by attaching the link (x, /x) to the path C xx at the 
end point x as is depicted in Fig. 1. Notice that Wa for a closed loop enters the r.h.s. of 
Eq. (|3.5| ) due to the presence of the delta-function. Therefore, the set of Eqs. (|3.3|) and 
).5\) is closed. 
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4 Properties of reduced and unreduced averages 



4.1 Reduced loop averages 

A prescription to calculate Wa(C) and G(C xy ) in the reduced theory can be obtained 
from Eq. fl2.6|) . The crucial role in this construction is played by the following integral 
over the unitary group: 

r b [$,^] = J dUe NtTWW ^%rt a Ut b tf , (4.1) 

where t a for a = 1, ... , N 2 — 1 are (hermitian) generators of SU(N) which are normalized 
by 

£ tr t a t b = 5 ab , and [t a ] kl = N5 a S kj - S^m (4.2) 

and I stands for the integral ( |2.14|) . Notice that the {/(l)-part is cancelled in this formula 
so that the integrals over U(N) and SU(N) coincide. 

It is convenient to consider I ab as a (iV 2 — 1) x (iV 2 — 1) matrix and 

$ a = itrt a $ (4.3) 

as a (iV 2 — l)-dimensional vector. Then, the counterparts of Wa and G in the reduced 
model read 

w Reduced = jLj^ ^jao (44) 

and 

QReduced = ^jL^ (4.5) 

where L is the length (in lattice units) of the appropriate contour. 

It is easy to see that W% educed and G Reduced given by Eqs. flO]) and (p]) satisfy for 
an arbitrary potential the same loop equation (|3.3|) as Wa(C) and G(C xy ) defined by 
Eqs. ( |1.3| ) and ( |1.4| ) in the original unreduced model. It is trivial to see, employing the 



representation flO ) in terms of the integral over U, that the equation coming from the 
shift Q3.2p of U coincides with Eq. ( |3.3D . The point is that we shift U only in one integral 
of the chain (what is an analog of the shift at one link) since one did not reduce U's (i.e. 
identify them). 

The situation is different, however, for Eq. (|3.5| ) which results from the shift ( |3.4|) of $. 
Now it is not possible to make any conclusion since the corresponding equation satisfied by 

Reduced 

can not be derived by a straightforward variation w.r.t. $ s . Fortunately Eq. (|3.5| ) 
can be exactly solved in the strong coupling region (see Section 5) and the result can be 
compared with the large mass expansion of G Reduced . 
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4.2 Large mass expansion 



The above conjecture about the equivalence of the reduced and original unreduced mod- 
els can be tested by comparing the loop averages at large values of the mass tUq. To 
calculate the large mass expansion, one needs the following expansion of I ab defined by 
the integral flPpfl: 

= § a m h + ... , (4.6) 

where . . . stands for the terms of higher powers in $ and \I/ which correspond to higher 
order of the large mass expansion. 

Applying this formula to the unreduced quantity (|1.4| ), one gets 
/ f d®(z)e- m o NtT ' s ' 2{z) ±tr$ 2 (z) \ 

«<c) - 11 ( J % (2) ,-*»JU» () ) + ■ ■ ■ = i^ L+ ' + ■ ■ ■ («) 

while the reduced prescription ( |4.5| ) gives 

QReduced = ^ tr $ 2 ^ L+1 + _ _ ^g) 

which coincides with Q4.7D since ■h tr Q 2 = Go- Therefore, one has explicitly demonstrated, 
in particular, the reduction to the leading order of the large mass expansion. 

A similar analyses of the adjoint Wilson loops leads to a different result. Let us now 
estimate the average on the r.h.s. of Eq. ( |1.3|) for a closed loop with L > 4 which does not 
contain parts that are passed back and forth. Exploiting the completeness condition (|4.2|) 
and Eq. ( f4.6|) , the result can be represented in the same form as ( |4.7|) but with an extra 
factor 1/N 2 : 

W A (C) = j^(G ) L + . . . . (4.9) 

This expression is proportional to 1/N 2 and vanishes in the large- N limit which is in 
agreement with general arguments of Ref . Q . The corresponding average in the reduced 
model can be calculated similarly to Eq. ( |4.8|) . The result is given by the r.h.s. of Eq. (|PD 
and now holds independently of whether the loop is closed or open. 

4.3 Intermediate coupling region 

When wig is decreased, the system undergoes phase transitions. According to the scenario 
of Ref. Q which is discussed in Section 1, a first order large- N phase transition should 
occur before the one associated with the continuum limit in order for the Kazakov-Migdal 
model to induce continuum QCD. While the strong and weak coupling phases always 
exist (the weak coupling phase is associated as usual with the perturbative expansion) 
the existence of such an intermediate phase has been only conjectured. Let us discuss 

2 The formulas for calculating the corresponding integrals over the Haar measure as N — > oo can be 
found in Ref. @. 
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consequences from the phenomenon of large-iV reduction in the intermediate coupling 
region. 

In the weak coupling region where the perturbative expansion is applicable the argu- 
ments of Section 2.2 about the large- N reduction work. 

The intermediate coupling region looks pretty similar to the weak coupling phase of 
Wilson lattice gauge theory. My arguments are based on the gauge field representa- 
tion fl2.18| ) of the Kazakov-Migdal model. Let us represent the induced action in the 
form 

S ind [U] = --J2/3 A (C)\tTU(C) | 2 (4.10) 
1 c 

where Pa(C) are some loop-dependent couplings which are determined by the potential 
V[$]. For the quadratic potential, one gets JI] 

MC) = (411) 

with /(C) being the length of the loop C. 

At N = oo the action Q4.10Q is equivalent to the following action 

S F [U] = -Nj2P(C)$ltrU{C) (4.12) 

c 

with the couplings /3(C) being determined by the self-consistency conditions 

P(C)=p A (C)W F (C;{P}) (4.13) 
where W F {C; {/?}) is the fundamental Wilson loop average 

W F (C ] {p}) = (±trU(C)} (4.14) 



and the averaging is taken with the action ( [4 .12 ). This procedure extends the one advo 



cated by Khokhlachev and the author |L8| for the single-plaquette adjoint action 

S A = -^ y £\trU(n)\ 2 (4.15) 



with the self-consistency condition given by Eq. (|4.13 ) with C = □ 



In the intermediate coupling region Eqs. (|4.13|) possess a nontrivial solution and the 
model can be described in terms of the fundamental induced action ( 4.12 ). On the other 



hand since the area law holds for the Wilson loop averages, I expect that the contribution 
of long loops to the induced action ( |4.12j ) is suppressed according to Eq. Q4.13| ). Therefore, 



properties of the intermediate coupling region should be similar to the standard lattice 
gauge theory. 



12 



5 An exact strong coupling solution 



5.1 Quadratic potential 

A drastic simplification of the loop equations (|3~3| ) and (|3.5f) is due to the fact f4| that 
in the strong coupling region (before the large-iV phase transition) the adjoint Wilson 
loops vanish in the large-iV limit except closed ones with vanishing area A m i n (C) of the 
minimal surface (i.e. contractable to a point owing to unitarity of t/'s): 

W A (C) = 5 QAmin{c) + O (jL} . (5.1) 

It can be shown that ( |5.1Q is consistent with loop equations. 

Given this behavior of adjoint Wilson loops entering the r.h.s. of Eq. ( |3.5| ), Eq. ( |3.5| ) 
admits a solution G(C xy ) = Gl depending on the single parameter L which is defined as 
the algebraic length of the loop (see ||IB|), i.e. the length after all possible contractions of 
the paths passing back and forth are made. For such an ansatz, Eq. ( |3.5|) can be written 

as 

2m\G L - G L _! - (2D - l)G L+1 = for L>1, 

2m 2 G - 2DG X = 1 . (5.2) 

Before solving this equation let me show of how it is satisfied in D — 1 where the 
exact result for Gl is given by the Laplace transform 



1 r oo \m - \fmt - 1 1 
- dae- m ° a I L (a) = ± . 7 (5.3) 

2 Jo 2Jml - 1 



with Il(ck) being modified Bessel functions. Eq. ( |5.2j ) can now be viewed, say, as the 
recurrent relation 

24(a) - I L _ x (a) - l L+1 (a) = (5.4) 

which is integrated according to Eq. (|5.3|) while 1 on the r.h.s. of the L = equation 
results from the fact that Io(0) = 1. 

Eq. (|5.2j ) can be solved by introducing the generating function 

oo 

G(X) = G L \ L (5.5) 

L=0 

with the result being expressed via the initial data, Go, by 

G(x) i{<G + D-l)X-(2D-l)G 

G{X) ~ 2mi\-V + l-2D ■ (5 ' 6) 
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While Go is defined as an average of -^tr$ 2 (:r) in the Kazakov-Migdal model, one 
can determine it from Eq. (|5.6|) by imposing the analytic structure of G(A) as a function 
of the spectral parameter A. For an arbitrary Gq the denominator in Eq. ( |5.6| ) has two 
roots and G(X) has therefore two poles. The exact solution at D = 1 given by Eq. ( |5.3| ) 
as well as the large mass expansion described in Section 4.2 lead to a G(X) possessing 
only one pole. This can be achieved in Eq. fl5.6|) by choosing Go in a proper way. This 
requirement unambiguously determines G(A) to be 



(2D 



G(X) = ; 

2 f m\(D - 1) + D^jml + 1 - 2D 

while Go is fixed to be 



ml + 1 - 2D - mi A + 2D - 1 



(5.7) 



Go = ° f . (5i 

m 2 (D - 1) + DJmQ + 1-2D 



I have checked that the solution ( |5.7| ) correctly reproduces the D — 1 solution flSlf ) as 
well as the leading order of the large mass expansion ([4.7|) at any .D. 

Some comments about the exact strong coupling solution (|5.7|) are now in order: 



Eq. ( |5.8|) precisely coincides with the result by Gross |J obtained by another method. 



While the solution ( |5.7| ) looks very simple, it corresponds in the language of the 
Kazakov-Migdal model to tedious calculations of integrals over unitary matrices 
with a subsequent averaging over $. In the reduced language ( |4.5| ) it remains still 
to calculate the integral ( |4.1| ) over unitary group and then substitute the value of 
the master field $ s . It would be very interesting to reproduce the result by this 
method. 

The solution ( |5.7|) is the first example of calculations of extended objects (since Gl 
is the average (|1.4j) for the loop of the length L in the lattice units) in the strong 
coupling region of the Kazakov-Migdal model while the results of Refs. |], 0, [5| refer 
to local objects like the spectral density. 

The very possibility to find such a simple exact solution to loop equations in the 
strong coupling region is related to a very simple form ( |5.1| ) of the Wilson loops. In 
particular, while the path C xy looks like a string in the strong coupling expansion 
of the standard lattice gauge theory, fluctuations of its shape are now suppressed. 
This is why the result depends only on L. 



The solution (|5.7p at D = 1 could be interesting from the viewpoint of matrix 
models of 2D gravity. 
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5.2 More general potentials 



While the solution ( |5.7| ) has been obtained for the quadratic potential, an analogous 
solution can be found for the more general potential 

iVtr~l/[$] = m 2 iVtr<l> 2 + #(tr$ 2 ) 2 (5.9) 

where g ~ 1 as iV — > oo to provide self-interaction of scalar field. This kind of potential 



was used |1(| in the matrix models in the context of D > 1 strings. 

The Kazakov-Migdal model with the quartic potential ( |5.9|) can be solved in the 
large- N limit for the following reasons. The model with the potential (|5.9|) is equivalent 
as N — > oo to the one with quadratic potential whose mass parameter rh 2 is defined by 
the self-consistency relation 

m 2 = ml + 2g G \^ (5.10) 

where Go I™, 2 * s gi ven by Eq. ( |5.8|) with replaced by rh 2 . Eq. (|5.10Q can be obtained 
naively replacing one tr $ 2 by the average value due to factorization. A rigorous proof 
of Eq. ( p.lOj) can be done using loop equations similarly to the proof [18] of the reduction 
of the adjoint action to the Wilson action which is discussed in Section 4.3. 

Eq. ( |5.10| ) can be used to study whether the large- iV phase transition occurs for the 
potential ( |5.9| ). Such a phase transition were occur if rh 2 would depend on Wq nonmono- 
tonically. A similar idea for the large- N phase transition to occur in lattice gauge theories 
was advocated in Ref. ]B5). For the scalar model of the type (|5.9p, it was employed in 



Ref. |T0] to obtain 7 string > 0. 



Given Eqs. (|5.10| ) and ( |5.8|) it is easy to calculate the derivative 

ftgj 1 _| 9(2D-l){D-l + ^^s) 
drfl2 (m 2 (D - 1) + DVrh* + l-2D) 2 

and see that it is positive for rh 2 > D where the gaussian model is stable. Therefore 
I conclude that, similarly to the case of the quadratic potential, there is no first order 
large-N phase transition for the Kazakov-Migdal model with the potential ( |5.9|) . 

An analogous study can be performed for the potential ( pL.5|) with an arbitrary function 
/ when the second term on the r.h.s. of Eq. ( |5.11p is multiplied by ~f"(Go). A conclusion 
is that one has to look for a more complicated potential. 
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Fig. 1 The graphic representation for <&{C xy ) (a) and $(C( x+ ^ x C xy ) (b) entering 
Eq. (|3.5|) . The bold points represent and $(x + /i). The (oriented) solid lines 
represent the path-ordered products U{C xy ) and U{C[ x+ ^ x C xy ). The color indices 
are contracted according to the arrows. 
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